Charge inhomogeneities due to smooth ripples in graphene sheets 
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We study the effect of the curved ripples observed in the free standing graphene samples on the 
electronic structure of the system. We model the ripples as smooth curved bumps and compute 
the Green's function of the Dirac fermions in the curved surface. Curved regions modify the Fermi 
velocity that becomes a function of the point on the graphene surface and induce energy dependent 
oscillations in the local density of states around the position of the bump. This effect is due to the 
change of the Pauli matrices with the position and is independent of the well known gauge potential 
generated by the spin connection. The corrections to the density of states are only due to the first 
effect and are estimated to be of a few percent of the flat density of states. Local probes such as 
scanning tunnel microscopy should be able to observe the predicted correlation of the morphology 
with the electronics. We discuss the connection of the present work with a recent observation of 
charge anisotropy in graphene and argue that the ripples can provide an alternative explanation. 
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I. INTRODUCTION 

Recent experiments of transmission electron (TEM) 
[l|, Q show that the suspended graphene samples ex- 
hibit an apparently random spontaneous curvature that 
can be visualized as ripples of various sizes that can 
reach a few angstroms high and several nanometers long. 
Observed ripples have also been reported later in scan- 
ning tunnelling microscopy I, H (STM). The issue of 
curvature of graphene and its possible influence on the 
electronic properties have been addressed before. In- 
spired by the physics of nanotubes and fuUerenes most 
of the works done on the electronic properties of curved 
graphene dealt with curvature induced by topological de- 
fects [1, H, H, Hi- In these works it was shown that 
conical singularities in the average flat graphene sheet 
induce characteristic charge anisotropies that could be 
seen in (STM) or (TEM). Charge anisotropies in mono- 
layer graphene have been recently observed in scanning 
single electron transistor experiments that show a 
distribution of electron-hole puddles at the Fermi level 
which could be responsible for the minimal conductiv- 
ity of graphene Charge anisotropies have also 
been reported in Electrostatic Force Microscopy exper- 
iments done in graphite [T^ . Charge inhomogeneities in 
graphene have been computed recently associated to im- 
purity states [l3| or elastic strain [ij] of a nature very 
similar to the one obtained previously with topological 
defects. It is by now clear that fluctuations in structure 
and charge are closely related and it seems that most 
types of disorder will induce charge anisotropies near the 
Dirac cone. 

In this work we study the effect of smooth curved por- 
tions of the type described in the experiments on the 
density of states of graphene by coupling the Dirac equa- 
tion describing the low energy electronic properties of 
graphene to a general metric describing a smooth curved 
piece in the average flat surface. The formalism is the 



same used in [8|, |9| to study the effect of pentagon and 
heptagon rings but here there are no curvature singu- 
larities. We obtain that the charge density couples to 
the shape of the graphene surface. The local density of 
states oscillates with the energy and the maximal correc- 
tion has a spacial extent of the order of the defect. The 
estimated value of the relative correction at the typical 
energies explored in the experiments can be of the order 
of a few percent of the flat value much bigger than the 
one obtained with elasticity models. The predicted cor- 
relation of the morphology of the sample with the elec- 
tronics should be observed in local probes as STM or 
TEM experiments. The present mechanism can explain 
qualitatively the recent observations of charge inhomo- 
geneities at zero energy found in graphene [lO| as due to 
the space variation of the Fermi velocity induced by the 
curvature. 

The article is organized as follows: Section ITTl explores 
the consequences of coupling the two dimensional Dirac 
Hamiltonian to a curved surface using a smooth gaus- 
sian bump as an example. At this level we see that as 
a consequence of the curvature two independent things 
happen: the Fermi velocity acquires a non trivial depen- 
dence on the position, and a gauge field is generated. 
Section IIIII is devoted to the computation of the local 
density of states of the system. We first describe the ef- 
fects of the curvature in terms of an effective potential 
and then compute the electron Green's function in the 
curved space to first order in perturbation theory. The 
perturbative parameter is related to the deviation of the 
graphene surface from the flat plane. Section IIVI ana- 
lyzes the results. In section|V]we discuss the possible ex- 
perimental consequences that can be extracted from this 
work in particular in connection with the anisotropies of 
ref. ■ In section IVII we present the conclusions and 
discussion. The appendix A is devoted to the definition 
and computation of the various geometrical factors asso- 
ciated to the problem. 



2 



II. TWO DIMENSIONAL DIRAC EQUATION IN 
A CURVED SPACE WITH POLAR SYMMETRY 

The most obvious way to study the effects of the topog- 
raphy of the sample on the electronic degrees of freedom 
is to couple the Dirac equation that governs the low en- 
ergy electronics to the curved surface. This ap pro ach has 
been applied to study curved fuUerenes [isl [T6I.IT7I. [isl. [T9| 
and to compute the response of electromagnetic charges 
to conical defects in planar graphene [20| in the frame- 
work of the equivalence between the theor y o f defects in 
solids and the three-dimensional gravity |2ll . [2^ . The 
coupling of the electronic degrees of freedom of planar 
graphene to conical defects has been explored in (sl. M. [23j . 
There it was found that a distribution of pentagons and 
heptagons induces characteristic inhomogeneities in the 
density of states of the graphene surface. In order to in- 
vestigate the effect of pure curvature on the electronic 
properties of graphene, the conical defects studied pre- 
viously present two difficulties. First they correspond to 
surfaces with zero intrinsic curvature; moreover the ex- 
trinsic curvature is accumulated at the apex of the cone 
where the surface has a singularity. It is then not clear 
if the results obtained - which look similar to the ones 
got with vacancies in [l3| - are due to the singularity or 
to the curvature. To disentangle the two effects we have 
studied the density of states of a flat graphene sheet with 
a smooth curved portion with intrinsic curvature. 

The massless Dirac equation in a curved spacetime is 
given by: 

i7^(5^ r!^)* = 0, (1) 

where = (p/^ ,vp^^), i=l,2. These curved 7 matrices 
satisfy the anticommutation relations 

{7^7''} = 25''''(x), (2) 

and in general become functions of the point in spacetime 
X = {t, r). ^^{x) is the spin connection of the spinor field 
that can be calculated using the tetrad formalism [23| and 
will be defined in the appendix A. 

What is needed is a metric describing the curvature 
of the surface. We will study the general case of a 
smooth protuberance fitting without singularities in the 
average flat graphene sheet. We start by embedding a 
two-dimensional surface with polar symmetry in three- 
dimensional space (described in cylindrical coordinates). 
The surface is defined by a function z{r) giving the height 
with respect to the flat surface z—Q, and parametrized 
by the polar coordinates of its projection onto the z=0 
plane. The metric for this surface is obtained as follows: 
We compute 

dz^ = (^fdr^ = af[r)dr\ (3) 

and substitute for the line element: 
ds^ = dr^ + r^dO^ + dz^ = (1 -I- af{r)) dr^ + r^dO^. (4) 




FIG. 1: (Color online) The smooth curved bump discussed in 
the text. 




FIG. 2: (Color online) Effect of the curved bump of fig. [T]on 
the local density of states of the graphene sheet. 



We will assume that the surface is asymptotically flat at 
long distances, so that / decays with r fast enough. We 
will also require / to go to zero sufficiently fast in r = 
so that the surface is smooth there. For clarity, we work 
out as an example the gaussian bump shown in Fig. [1] 
defined by: 

z = Acxp(-rVfc^), (5) 

so that 

dz^ = -^4r2 exp{~2ryb^)dr\ (6) 

which corresponds to eq. ^ with 

a = {A/bf , /(r) = 4(r/6)2exp(-2rVfo^). 

The ratio of the height to the mean width of the gaus- 
sian will be our perturbative parameter. Values of b of 
the order of 0.1 to 0.3 times the height A give bumps 
of magnitude and shape comparable with the ripples re- 
ported in [U,!!!]. 

Since we are going to work with the Dirac equation, 
our calculation is formally relativistic and the full space- 
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time metric must be used. The line element reads: 

ds^ =dt^ - {l + af)dr^ ~r^d0^, (7) 
from where we can write the metric in a more usual form: 



10 

.g^.- I -(l + a/(r)) 



(8) 











Since we are dealing with a problem with axial symmetry, 
we will work in polar coordinates. The Dirac Hamilto- 
nian of the plane (flat) in polar coordinates can be writ- 



J 



ten as 



Hflat = ^Vp 



I r 2r 



1 



2r 



(9) 



where, as discussed in Appendix A, a constant spin 
connection has been generated that is a "pure gauge" 
and can be rotated away by a different choice of local 
coordinates [3^. The calculation of the different geomet- 
ric factors for the metric ([7]) is given in the appendix A. 
Adding them together we get the curved Hamiltonian 



TiVf. 







(10) 



r 



where the effective gauge potential is related to the coef- 
ficient of the spin connection computed in the appendix 
(IMD) by 



l-(l + a/)-^/^ 

Aa — = . 

2r 2r 



(11) 



Comparing ([9]) with (fTO|) we can read that the curved 
bump has produced an effective Fermi velocity Vr in the 
radial direction given by 



Vr{r,e)=VF[l + af{r)r^''', 



and an effective magnetic field perpendicular to the 
graphene sheet given by 



1 



drirAg 



1 



4r (1 a/)3/2 ■ 



(13) 



The magnitude of this effective magnetic field is esti- 
mated to be of the order of 0.5 to 2-3 Tesla in the region 
spanned by the bump, compatible with the estimations 
given in [25], and it will play the same role in the is- 
sue of the weak localization of graphene as the effective 
magnetic fields discussed there and in 

We note that in general the effective Fermi velocity 
will be smaller in magnitude than the free one. For a 
general curved surface described in polar coordinates by 
z = z(r), the effective Fermi velocity will be 



Vr 



(14) 



In a most general case we will have the two components 
of the velocity changed but always to a smaller value. 

In the next section we specify the method to compute 
the density of states through the electron Green's func- 
tion. For this purpose it is more convenient to follow the 



appendix and rewrite the Dirac equation in the form 

^- = 0, (15) 



i7°5o + iT{e)dr + tT'{0)— + Vir, ( 

r 



which is the flat Dirac equation in polar coordinates in a 
sort of potential V given by: 



V{r,0)=iT{e) l-(l + a/)-^/^ iir-^r). (16) 

2r 



1 



(12) where r{0) = 7^ cos 6 + 7^ sinfl. This effective potential 
will be used in the next section to compute the local 
density of states of the system. We can read in it two 
different terms related with our previous discussion: the 
derivative term has its origin in the effective r-dependent 
Fermi velocity and the term proportional to 1/r comes 
from the effective gauge field. 



III. THE GREEN'S FUNCTION IN A CURVED 
SPACETIME. APPROXIMATIONS 

We drop the polar coordinates for a moment to outline 
the procedure to obtain the Dirac propagator in a sort of 
perturbative expansion. 

Since we will be interested in small deformations from 
the flat membrane, we will compute the first order cor- 
rection to the "flat" propagator in the small parameter 
a. In our example, a = {A/b)"^ is the (squared) height 
to length ratio of the gaussian, so for typical ripples in 
raphene a « 0.01, since this ratio is of the order of 0.1 

The equation for the exact propagator in the curved 
space-time is: 
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{d^ + n,,) G{x,x') = S{x-x'){-gy 



(17) 
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where is the spin connection computed in the ap- 
pendix (|A21[) and \/—g is the determinant of the metric 
given in (jA7|l . Since the flat propagator equation is re- 
covered when a = 0, expanding the left hand side to first 
order in a will give this flat equation plus a first order 
general term that we will call V. We can as well expand 
(— g)"^/^ = 1 — af{x), and sending the f term to the left 
hand side we get an equation resembling the flat prop- 
agator equation in a sort of potential generated by the 
metric. 

(^7^9^ + V)G{x, x') + af{x)S{x - x') = S{x - x'). (18) 



Next we need the Dirac propagator in polar coordi- 
nates. We will get it by noting that the Dirac and Klein- 
Gordon propagators are related by: 



G% — D{x)G%Q, 



(22) 



where D(x) is the Dirac operator. The Dirac operator in 
polar coordinates is: 



T{e)dr + T'{e)— 

r 



(23) 



This equation can be solved by the usual pertubative 
expansion of G in a potential. Note that two approxima- 
tions to order a are taking place: first, we expand the 
curved space exact equation for G to make it resemble 
the flat equation in a potential, and then we use a per- 
turbative expansion in the potential to get the first order 
correction to G. 

The solution for G is given by: 

Gix, x') = Go{x - x') - Ga{x - x')af{x) + 

dx"Go{x, x")V{x")Go{x", x'). (19) 



where we have defined: 



r(e') 7^ cos(6l) + 7^ sm{9) , 
r'{e) = -7^ sin(6l) + 7^ cos(6'). 



(24) 



The Klein-Gordon propagator in polar coordinates is 



GKG{r'r\E)^—HoiE\r-r'\), (25) 

where Hq is the zeroth order Hankel function ^J\. Ap- 
plying eqs. ([23|) and ([22]) we get: 



We now proceed to use this expansion in our particular 
case. The determinant of the metric is just: 



_ il + c^f)-'/' 



^ 2 



Expanding (jl6p to first order in a we get: 



V{r,0) = ir{9) 



(20) 



(21) 



Noting that the flat S function in polar coordinates is 
5{r — r')/r, we can use eq. in polar coordinates to 
get the first order correction to the propagator. 



GD{T^r\E) = ^^'^H^{E\ 
E 



^—^^H,{E\r-r'\)[rT{9)-r'T{e')]. 



(26) 



Now we just use eqs. ((T^ and (f^^ . and the first order 
potential ((2T|) to get the local density of states 



p{E,v) ^ ImTr[G{E,Y,Y)-i\ 



(27) 



After two partial integrations we obtain the local density 
of states as: 



p{r',E) = g 



E 



27r [vpfif 



af{r') 



a 



drde. 



r ~ r' cc^ H — 0' 
[r^ + r'-^ — 2rr' cosiO — O'))^''^ 



(28) 



where g is the spin and valley degeneracy, A(r) — \i — r' 



and Yi , are Bessel functions 
result in the next section. 



27| . We will analyze this 



IV. RESULTS 



The main result of this work is the confirmation of 
the fact that the morphology of the graphene samples is 
correlated with the electronic properties. In particular 
the presence of ripples, irrespective of their origin, in- 
duces corrections to the density of states and will affect 
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FIG. 3: (Color online) Correction to the density of states 
(dotted line) in arbitrary units, versus the shape of defect 
(red - upper - line) and curvature of the defect (yellow - 
middle- line) for a gaussian bump of an average width of 50 

A. 



the transport properties. The analysis of this work is 
quite general and the results will apply to any smooth 
surface. We have considered as an example a gaussian 
bump of an average spacial extent of 2 to 20 nanome- 
ters similar to the ripples described in [l], [2^ . As can be 
read from eq. (|28|) due to the Bessel functions the correc- 
tions to the local density of states (LDOS) show spacial 
oscillations whose frequency grows with the energy and 
whose amplitude decays with the distance as (1/r). Fig. 
[5]shows the correction to the LDOS induced by the shape 
of Fig. [T]with a mean width b of 50 A and for an energy 
of i? = 0.1 eV. The color scale is indicated in the fig- 
ure. Lighter (darker) color means a positive (negative) 
contribution with respect to the flat graphene sheet at 
the given energy. The maximal value of the correction 
related to the bare LDOS p{E, r)/po{E,r) for a bump of 
a ratio A/b ~ 0.1 is of the order of a 1 percent. If the 
height of the ripples goes up to 0.3, the maximal value 
of the correction to the LDOS due to the curved portion 
can reach a 10 percent. Larger values are possible but 
they would correspond to higher values of our perturba- 
tive parameter a and made the result questionable. The 
bare LDOS being proportional to the energy (E), allows 
to relate the correction obtained with the energy at which 
the experiment takes place. 

At low energies the maximal absolute value of the cor- 
rection is correlated with the zero of the curvature which, 
in the particular case of the gaussian bump, coincides 
with the mean width \^b. Fig. [3] shows the correction 
to the local density of states (dotted line) in arbitrary 
units, versus the shape of defect (red - upper ~ line) and 
the curvature of the defect (yellow -middle- line) for a 
gaussian bump of an average width of 50 A. The energy 
is 0.1 eV. The figure represents real space in polar coor- 
dinates. The horizontal axis is the r coordinate while the 
vertical axis represents real height in the case of the up- 
per line (shape of the gaussian bump considered). The 
middle line (curvature) given in eq. (IA14P is not mea- 
sured in units of length. Finally the lowest curve gives 
the correction to the LDOS in arbitrary units. This re- 
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FIG. 4: Plot of the location of the maximal correction of the 
LDOS for a fixed energy as a function of the extension of the 
Gaussian b. 



suits seems to be at odd with related works based on 
topological defects [H, [2^ or elasticity [13] that corre- 
late the physical effect of curved portions with the actual 
value of the geometrical curvature. 

We notice here that of the two effects of the curvature 
discussed in this work, only the effective magnetic field 
coming from the spin connection can be compared with 
previous works. In our case and in the general situation 
of having a smooth shape with axial symmetry, the effect 
of the effective gauge field vanishes at first order in per- 
turbation theory and the corrections to the local density 
of states come exclusively from the spacial dependence 
of the Fermi velocity. This effect has not been noticed 
before because most of the previous works coupling the 
Dirac equation to curved space dealt either with spheri- 
cal shapes where the correction to the Fermi velocity is 
constant and can be scaled out, or with conical shapes 
whose intrinsic curvature is accumulated at the apex sin- 
gularity. 

A discussion on the various gauge fields that arise in 
the physics of graphene and their physical consequences 
will be published elsewhere p8| . 

For simplicity we have modelled the ripples with shapes 
that are axially symmetric and the axial symmetry is 
explicit in the results. More general shapes would made 
the calculation much more complicated without altering 
the main results. 

We can observe in figs. [2] and [3] that the maximal value 
of the correction to the LDOS is concentrated in the re- 
gion spanned by the bump. In an STM experiment done 
on the curved surface the results on the LDOS plotted 
on the flat surface should resemble the rings of fig. 
and would be in a very precise correspondence with the 
morphology of the sample. 

To verify the apparent pining of the maximal value of 
the correction to the zero of the curvature observed in [3] 
we have explored a set of gaussian shapes with different 
widths. The result is plotted in fig. 21 The correlation 
is good for the physical values of the bump from 2 to 
5 nanometers (50 A). Fig. [5] shows the dependence of 
the correction to the LDOS with the energy. The four 
plots show the correction induced by the same bump of a 
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FIG. 5: Correction to the density of states of a bump of fixed 
width (50 A) for different values of the energy. 



width of 5 nanometers for growing values of the energy. 



V. POSSIBLE EXPERIMENTAL 
CONSEQUENCES. 

We will here discuss the possible connection of this 
work with the recent experimental observations of 
electron- hole puddles of ref . [l^ [sH] . 

The main experimental facts that we have to address 
are the intrinsic disorder length scale of approximately 
30 nm and the magnitude of the density fluctuations of 
an average value An ~ ±4.10^°c77i~^ that can reach a 
maximum value of lO^^cm^^. 

In the experiment described in p^ . the local inverse 
compressibility d^/dn is measured as a function of a 
back-gate voltage. In the case of having a flat graphene 
sheet this quantity can be computed analytically using 
the free dispersion relation and density of states, and it 
becomes more complicated when the latter is corrected 
due to the curvature. 

Before entering in a more detailed a naly sis we can 
explain qualitatively the observation of as follows: 
From our result of Section IIIII we note that a variable 
correction to the density of states of the type shown in 
causes a global, constant variation of the Dirac point 
Eu with respect to the charge neutrality point, because 
of the redistribution of the electronic charge when the 
density of states is modified locally (notice from Fig. [3] 
that the corrected LDOS is always lower than the flat 
one). This mechanism is the one proposed in ref. 14] in 
a different context. An observation like the one done in 
[l0| will see electron and hole puddles of the size of the 
ripples without the need to argue for charged impurities 
or other kinds of disorder. To give a number, we can 
estimate the density of carriers that we would need to 
induce to reach the Dirac point at a given position of the 
sample in this case as the integral of the correction of the 



density of states from some lower cutoff (of the order of 
-0.7 eV ) to 0. In the case of the maximal correction we 
find this density to be of order An ~ ilO^^cm"^. 

For a more quantitative discussion we consider a flat 
graphene sheet with some general Eo measured with re- 
spect to the charge neutrality point Ecnp, which we set 
to zero. The local density of states would be given in 
such a case by: 



p{E,r) 



\E-ED{r) 



9s 



2tt 



[flVp) 



(29) 



where gs—"^ is the degeneracy (spin and valley). When 
a voltage is applied, a density of carriers is induced lo- 
cally, and thus the local chemical potential varies as well, 
though not linearly with the voltage. The induced num- 
ber of carriers in this case can be computed as: 



dE 



\E-Er 



9S 



2ti (hvpy 



dE 



\E\ 



9s 



27r {hvpY ' 
(30) 

The integration yields, taking appropriate care of the 
signs due to the absolute value: 

[sign{-ED + m)(--Ed + ^if - sign{-ED)E'U gs 

fl — — . 

(31) 

The rate of change dn/dfj, is again p{—E£i + /i), and 
solving for /i in (I3ip . substituting the value and inverting, 
we obtain: 



dp 
dn 



hvf 



ngs 



i-rr'-' {hv f)'-' 



1/2 ' 



(32) 



which reduces to the formula given in ref. [10| when 
Ed = 0. We can see that if Ed is shifted from zero, the 
curve for the inverse compressibility keeps its shape but 
it is displaced in the n axis, as shown in the experiment. 

However there is an alternative interpretation of the 
displacement of the curves. Instead of having a variable 
Ecnp from point to point, we could explain the shift by 
having a variable as predicted in this work. We con- 
sider now the curved case, in which in a qualitative first 
approximation we just use eq. (|32p where Vp vpir). 
We can see in eq. ((32|) that this substitution has two ef- 
fects: first, it makes dp/dn to scale with vp{r), and also 
the divergence is displaced by an amount proportional 
to {EolvF{r)Y , different from the flat case. The scaling 
effect provides an experimental test of the variable Fermi 
velocity. Extracting the value of VF{r) from the shape of 
the curves at each point we can get a corrected map for 
Ed- We have analyzed carefully the experimental data 
of r ef. [13] m to check this issue but the precision is not 
enough to deduce the change in the slope of the curves 
from point to point. Perhaps this effect can be tested in 
angle resolved photoemission experiments. 
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A related comment concerns the spatial variation of 
the Fermi velocity obtained in this work. Different rip- 
ple sizes randomly distributed over the sample will give 
rise to a landscape of Fermi velocities that will directly 
account for the size of the inhomogeneities reported. We 
note that the effective Fermi velocity is the (only) fitting 
parameter used in [To| . 



VI. CONCLUSIONS AND DISCUSSION 

We have studied the effect of curvature on the elec- 
tronic properties of a graphenc sheet as due to the spinor 
nature of the electron wave function. With complete gen- 
erality we have seen that the curvature of the graphene 
sheet has two distinct effects with physical consequences. 
First the covariant derivative induces an effective mag- 
netic field (the spin connection) that depends on the 
shape of the curved surface. This effect has been dis- 
cussed before in the literature and can be understood in 
many ways. The second effect is even more interesting 
and comes from the curved Pauli matrices. They modify 
the Fermi velocity that becomes a function of the point 
on the graphene surface and is always lower than the free 
velocity. This is not to be confused with the renormaliza- 
tion of the Fermi velocity induced by electron-electron in- 
teractions 29] or disorder [3, [s^l where the Fermi velocity 
becomes effectively a function of the energy. This effect 
has not been noticed before because most of the previous 
works coupling the Dirac equation to curved space dealt 
either with spherical shapes where the correction to the 
Fermi velocity is constant and can be scaled out, or with 
conical shapes whose intrinsic curvature is accumulated 
at the apex singularity. 

To model the ripples observed in graphene we have ap- 
plied the formalism to study a smooth surface with local 
curvature that is asymptotically flat. We have studied 
the changes in the density of states by computing the 
electronic Green's function in the presence of the curved 
bump. The effective potential generated by the curvature 
induces oscillations in the local density of states that af- 
fect significantly the region spanned by the bump. The 
potential couples to the electronic Dirac equation as a 
gauge field with two parts. One can be identified with an 
effective non abelian magnetic field of the type described 
previously associated to topological defects or external 
charges. This part does not contribute to the LDOS to 
first order in or perturbative expansion due to the traces 
over the gamma matrices. It will nevertheless affect other 
observables and will generate a density of states in non- 
perturbative calculations of the type recently presented 
in [3l|. The second piece of the effective potential repre- 
sents a gauge field with a derivative coupling and comes 
from the modified Fermi velocity. This is the one pro- 
ducing the results presented in this work. A conservative 
estimate gives a relative correction to the fiat density of 
the order of 0.5 to a few per cent at the energies usually 
explored with STM probes that should be able to corre- 



late the morphology with the predicted correction. The 
estimated corrections due to curvature are considerable 
larger than the ones recently obtained with the theory of 
elasticity fl^. They are also different: in most cases due 
to the traces over the gamma matrices, the correction 
to the density of states that would come from the spin 
connection vanishes and what remains arises exclusively 
from the derivative term in the potential (the correction 
to the Fermi velocity). The space variation of the Fermi 
velocity induced by the curvature provides an alterna- 
tive explanation for the charge anisotropics observed in 
and reinforces the effect. We propose a mechanism 
to disentangle the two effects (a space dependent chem- 
ical potential versus a space dependent Fermi velocity) 
that would provide an experimental test of the curvature 
effects described in this work. 

The mechanism proposed based on the variable Fermi 
velocity due to the curvature of the sheets can be distin- 
guished from other mechanisms by the fact that a vari- 
able vp in addition to displace the position of the curves 
would change their slope at different positions of the sam- 
ple providing an experimental test of the model discussed 
in this work. 

From the discussion done in the paper is clear that 
a part of the described modification of the Fermi veloc- 
ity can be transferred by a gauge transformation to the 
vector potential and viceversa. Any formulation will pro- 
duce the same results when computing observable quan- 
tities. In this respect we should not identify the effective 
Fermi velocity with any tight binding parameter. As hap- 
pens in general relativity, curvature means interactions. 
The Fermi velocity measured in the photoemission ex- 
periments or by any other means is the result of the bare 
value and the interactions so there is no paradox here. 

The presence of a random distribution of curved por- 
tions as the one discussed in this work will affect the 
transport properties of the system. In this case the ef- 
fective magnetic field will p lay a role similar to the one 
discussed in the literature This issue is currently 

under investigation and it will be reported in a different 
work. 
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APPENDIX A: THE GAUSSIAN BUMP. 
GEOMETRIC FACTORS. 

The behavior of spinors in curved spaces is more com- 
phcated than that of scalar or vector fields because their 
Lorentz transformation rules do not generalize easily to 
arbitrary coordinate systems. Instead of the usual metric 
ff/ij/ we must introduce at each point X described in arbi- 
trary coordinates, a set of locally inertial coordinates 
and the fielbein fields e° (a;), a set of orthonormal vectors 
labelled by a that fixes the transformation between the 
local and the general coordinates: 



el[X) 



(Al) 



We will later compute the fielbein for our particular met- 
ric. The curved space gamma matrices 7^(2;) satisfying 
the commutation relations 



Despite the fact that the spin connection appears to be 
non trivial, the Riemann curvature which is the "observ- 
able" quantity and does not depend on the choice of co- 
ordinates is zero as corresponds to fiat space. 
The fielbein fields e"^ satisfy: 



(AlO) 



where rjab is the identity matrix in two dimensions. 

This relation does not fix e"^ uniquely. There are two 
natural choices: one is 



1 
r 



and the other one is 



COS0 — rsin6' 
sin 9 r cos 9 



(All) 



(A12) 



{7^7j = 2g^„ (A2) 
are related with the constant, flat space matrices 7° by 
^^{x) = e^^7a. (A3) 
The spin connection fi^ (x) is defined from the fielbein by 
1 



with 



T^A a 



(A4) 



(A5) 



The two choices can be visualized as associated to fiat 
local frames that at points of constant r have fixed di- 
rections (last) or rotate with the polar angle (first). The 
first choice leaves the gamma matrices as in the cartesian 
plane and induces a constant gauge connection whose 
"associated magnetic field" is obviously zero. 

The second choice transforms the fiat gamma matrices 
and does not induce a gauge connection. 

Let us now compute the geometric factors related with 
the metric of eq. ([8|): The affine connection F^^ for the 
metric (151) is: 



2(1 + a/)' 



r yO 
l+af' ^rl 



,(A13) 



where T^^ is the usual affine connection which is related 
to the metric tensor by 33] 



1_ v\i^9av_ 

2^ ^ dxf^ 



dx° 



dx" 



}■ (A6) 



Finally, the determinant of the metric needed to define 
a scalar density lagrangian is given by 



^= [det(5^^)] 



1/2 _ 



detKix)]. 



(A7) 



Before going to the computation of the geometric factors 
related with the metric of eq. ([5]) we will apply the for- 
malism to the fiat space in polar coordinates what will 
help to clarify the physical discussion later. The two di- 
mensional metric of the flat space in polar coordinates 
is 



1 

r2 



(A8) 



where /' = df /dr, and the rest of the elements are zero 
or related by symmetry. 

The geometrical (gaussian) curvature K of the shape 
given by eq. ^ is 



af{r) 



2r{l + af{r)y 



The fielbein fields e"^ satisfy: 

where g^i, is our metric given in eq. ([5]) and 

1 
?7a6 = I -1 
0-1 



We choose the e" to be 



(A14) 



(A15) 



(A16) 



The affine connection F^^ that only depends on the met- 
ric is 



1 — U, i — ' , L re 



1 



(A9) 



e\ = e°e = 

eV = (1 + a/)^/^cos6l e^^ = -rsin( 
e\ = {l + afY/^sm9 e^g^rcos9- 



(A17) 
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that reduce to the flat set (|A12p when a = 0. Now we 
can compute the spin connection cocflicients, 



which are found to be: 



(A18) 



(A19) 



the rest being zero or related by symmetry (the s pin con- 
nection u) is antisymmetric in the upper indices [33|). 
The spin connection 



(A20) 



turns out to be 

fij = 0, ?lr = 0, ng 



Finally the Dirac equation coupled to the curved surface 
is: 



(A22) 



Substituting all previously computed elements and with 
some more algebra, we can cast ([T} into the form: 



ij°do + iT(e)dr + ir'(e)— + v{r, i 

r 



* = 0, (A23) 



which is the flat Dirac equation in a sort of potential V 
given by: 



1 
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